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ABSTRACT 

AT / ^ 

The three dimensional three degree of freedom re-entry equations of 
motion and the input requirements for their solution on the IBM 7090 
digital computer are derived. The material presented is the first formu- 
lations required to achieve eventually a hardware error analysis capability 
for three dimensional six degree of freedom closed loop re-entry guidance 
and control. The presentation consists of the derivation of the three 
dimensional three degree of freedom re-entry trajectory expressed in terms 
of terrestrial navigational coordinates. The derivation of the coordinate 
transformation equations required to express the solution of the above 
equations in terms of inertial platform coordinates is also presented. 
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DEFINITION OF SYMBOLS 


SYMBOL 

DEFINITION 

UNITS 

A I 

Space-fixed azimuth of initial space-fixed velocity 

rad 

\ 

Earth-referenced azimuth of earth-referenced wind 
velocity 

rad 

A z 

Earth-referenced azimuth of earth-referenced aero- 
dynamic velocity 

rad 

C D 

Aerodynamic drag coefficient for zero angle of attack 

rad 

o 

C D 

Total aerodynamic drag coefficient 

rad 

C L 

Total aerodynamic lift coefficient 

rad 

D 

Magnitude of the total aerodynamic drag force 

kg 

F a 

Total aerodynamic force 

kg 

? D 

Total aerodynamic drag force 

kg 

? G 

Total gravitational force 

kg 

? L 

Total aerodynamic lift force 

kg 

G 

Universal gravitational constant 

m 4 /kg-s 

G 

r 

Magnitude of the gravitational acceleration in the 
geocentric radial direction 

m/ s 2 


Magnitude of the gravitational acceleration in the 
direction of increasing earth longitude 

m/s 2 

% 

c 

Magnitude of the gravitational acceleration in the 
direction of increasing geocentric latitude 

m/s 2 

J 

Second order figure constant of the earth 

rad 

L 

Magnitude of the total aerodynamic lift force 

kg 

M 

Mass of the earth 

kg-s 2 /m 

R 

Displacement from the center of the earth 

m 

S 

Vehicle effective aerodynamic surface area 

iv 

m 2 


DEFINITION OF SYMBOLS (Cont 'd) 

SYMBOL DEFINITION UNITS 

U Gravitational potential of the earth m 2 /s 2 

Aerodynamic velocity m/s 

V A Magnitude of aerodynamic velocity m/s 

Magnitude of initial space-fixed velocity m/s 

Magnitude of earth -referenced wind velocity m/s 

a^ Empirical aerodynamic drag coefficient caused by rad 

vehicle angle of attack 

a^ Empirical aerodynamic lift coefficient caused by rad 

vehicle angle of attack 

a Magnitude of space-fixed acceleration in m/s 2 

instantaneous direction of increasing geocentric 
displacement 

a^ Magnitude of space-fixed acceleration in instantaneous m/s 2 

direction of increasing space-fixed longitude 

a^ Magnitude of space-fixed acceleration in instantaneous m/s 2 

c direction of increasing geocentric latitude 

b^ Empirical aerodynamic drag coefficient as a function rad 

of angle of attack squared 

b Empirical aerodynamic lift coefficient as a function rad 

of angle of attack squared 

h Q Initiai geocentric altitude m 

h Geocentric altitude m 

j Inclination of the initial orbit plane rad 

m Vehicle mass kg-s 2 /m 

q Dynamic pressure kg/m 2 

r Magnitude of initial radial displacement from the m 

° center of the earth 


v 



DEFINITION OF SYMBOLS (Cont'd) 


SYMBOL 

DEFINITION 

UNITS 

r 

Magnitude of radial displacement from the center of 
the earth 

m 

r E 

Equatorial earth radius 

m 

r 

P 

Polar earth radius 

m 

• 

r 

0 

Magnitude of initial radial velocity from the center 
of the earth 

m/s 

• 

r 

Magnitude of radial velocity from the center of the 
earth 

m / s 

• • 
r 

Magnitude of radial acceleration from the center of the 
earth 

m/ s 2 

t 

Time 

s 

V 

a E 

Space-fixed velocity of air molecule in a stationary 
atmosphere 

m/s 

V 

a 

• w 

Space-fixed velocity of air molecule in flowing 
atmosphere 

m/s 

V 

r 

Magnitude of space-fixed velocity in the instantaneous 
direction of increasing geocentric displacement 

m/s 

V A 

Magnitude of space-fixed velocity in the instantaneous 
direction of increasing space-fixed longitude 

m/ s 

% 

c 

Magnitude of space-fixed velocity in the instantaneous 
direction of increasing geocentric latitude 

m/ s 

X 

Reference space-fixed earth-centered cartesian 
coordinate in equatorial plane 

m 

X 1 

Earth-centered cartesian coordinate generated by first 
Euler rotation 

m 

X 

2 

Earth-centered cartesian coordinate generated by 
second Euler rotation 

m 

y 

Reference space-fixed earth-centered cartesian 
coordinate in equatorial plane 

m 


vi 




DEFINITION OF SYMBOLS (Cont'd) 


SYMBOLS 

DEFINITION 

UNITS 

yi 

Earth-centered cartesian coordinate generated by first 
Euler rotation 

m 

y 2 

Earth- centered cartesian coordinate generated by 
second Euler rotation 

m 

z 

Reference space-fixed earth -centered cartesian 
coordinate in equatorial plane 

m 

Z 1 

Earth-centered cartesian coordinate generated by first 
Euler rotation 

m 

z 2 

Earth-centered cartesian coordinate generated by 
second Euler rotation 

m 

V 

"Del" vector operator, or "Nabla 11 

l/m 

At 

Numerical integration time step 

s 

(At) p 

Data printout time step interval 

s 

A 

S pa c e - f ixed long i t ude 

rad 

A 

Space-fixed longitudinal angular velocity 

rad/s 

A 

Space-fixed longitudinal angular acceleration 

rad/s 

0 

Space-fixed total angular velocity 

rad/s 

a 

Angle of attack 

rad 

7 

Bank angle: defined as the angle between the 
aerodynamic lift vector and the vertical geocentric 
plane containing the aerodynamic velocity vector 

rad 

5 

Space-fixed displacement in the initial cross range 
direction 

m 

T] 

Space-fixed displacement in the direction of the 
initial geocentric altitude 

m 


Displacement of sea level surface from the center of 
the earth in the direction of the initial geocentric 
altitude 

m 


vii 



DEFINITION OF SYMBOLS (Cont'd) 


SYMBOL 

DEFINITION 

UNITS 

\ 

Displacement from sea level surface in the direction 
of the initial geocentric altitude 

m 

0 A 

Earth-referenced dive angle of aerodynamic velocity 

rad 

6 I 

Space-fixed dive angle of initial space-fixed velocity 

rad 

X 

o 

Initial earth longitude 

rad 

X 

Earth longitude 

rad 

X 

o 

Initial earth longitudinal angular velocity 

rad/ s 

X 

Earth longitudinal angular velocity 

rad/s 

• • 

X 

Earth longitudinal angular acceleration 

rad/s 2 


Space-fixed displacement in initial downrange direction 

m 

p 

Atmospheric mass density 

kg-s 2 /n 

cl) 

c 

Initial geocentric latitude 

rad 

o 

<fr 

c 

Geocentric latitude 

rad 

* 

c 

0 

Initial geocentric latitudinal angular velocity 

rad/ s 

i 

c 

Geocentric latitudinal angular velocity 

rad/s 

V 

c 

Geocentric latitudinal angular acceleration 

rad/s 2 

(j 0 

Angular velocity of earth rotation 

rad/s 
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A DERIVATION OF THE THREE DIMENSIONAL THREE DEGREE OF FREEDOM 
RE-ENTRY EQUATIONS OF MOTION IN TERRESTRIAL NAVIGATIONAL 
COORDINATES WITH SUBSEQUENT RESOLUTION INTO INERTIAL 
PLATFORM COORDINATES 

By Carl L. Colwell 
SUMMARY 

The three dimensional three degree of freedom re-entry equations of 
motion and input requirements for their solution on the IBM 7090 digital 
computer are derived. The material presented is the first formulations 
required to achieve eventually a hardware error analysis capability for 
three dimensional six degree of freedom closed loop re-entry guidance and 
control. The presentation consists of the derivation of the three simulta- 
neous differential equations of motion required to solve for the three di- 
mensional three degree of freedom re-entry trajectory expressed in terms 
of terrestrial navigational coordinates. The derivation of the coordinate 
transformation equations required to express the solution of the above 
equations in terms of inertial platform coordinates is also presented. 

SECTION I . INTRODUCTION 

The problem of hardware error analysis requires that accurate simula- 
tion models of various types of mission trajectories be available. These 
trajectories are used to determine the effects of various instrumental 
errors on the overall accuracy of the guidance systems. 

The general approach to obtaining three degree of freedom equations 
of motion for a three-dimensional re-entry trajectory is a treatment of 
classical particle dynamics utilizing vector kinematics for the derivation 
of the space-fixed acceleration. The dynamical model of the re-entry 
vehicle is a point mass subjected to the forces of the gravitational field 
of the earth and aerodynamic lift and drag. 

First, the space-fixed acceleration, expressed in terms of a moving 
coordinate system, is derived by vector kinematics for direct substitution 
into Newton's Second Law. The gravitational and aerodynamic forces are 
then expressed in terms of that moving coordinate system and equated to 
the vehicle mass times the space-fixed acceleration. Subsequent cancella- 
tion of the unit vectors from the three component equations then yields 
the three scalar differential equations to be solved simultaneously for 
the three dimensional re-entry trajectory. 
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A judicious sequence of Euler angle rotations of the moving coordinate 
axes then transforms the trajectory parameters expressed in the moving 
coordinate system into equivalent parameters in a coordinate system 
oriented in the direction of an idealized inertial platform coordinate 
system. 


SECTION II. DEFINITION OF COORDINATE SYSTEMS 


Two rectangular coordinate systems and one spherical coordinate 
system are utilized in the treatment of the re-entry problem for which 
the center of the rotating earth is considered as being space-fixed. 

These three coordinate systems are illustrated in Figure 1. 

The l x , ly, l z unit vectors represent a right-handed space-fixed 
rectangular coordinate system which has its origin at the center of the 
earth. The x and y axes are oriented in the equatorial plane with the 
x axis oriented in a reference space-fixed meridian plane, and the z 
axis collinear with the North Pole. This coordinate system serves only 
as a visual aid in deriving the space-fixed acceleration and is not 
intended to convey any other physical significance. 

The , 1^, 1^ unit vectors represent a rotating spherical coordinate 
c 

system which also has its origin at the center of the earth. The unit 
vectors are oriented in the right-handed sense and represent the 
instantaneous direction of increasing geocentric latitude, the instan- 
taneous direction of increasing geocentric displacement, and the 
instantaneous direction of increasing space-fixed longitude, respectively. 
The three simultaneous differential equations of motion are derived in 
terms of this coordinate system. 


The 1^, 1 , and 1^ unit vectors represent a space-fixed coordinate 

system which has its origin on the surface of the earth directly beneath 
the initial re-entry point. The unit vectors are oriented in the right- 
handed sense and represent the initial downrange direction, initial 
geocentric altitude direction, and initial cross range direction, respec- 
tively. The solution of the equations of motion is resolved into this 
coordinate system. 


SECTION III. DERIVATION OF SPACE-FIXED ACCELERATIONS 

To derive the space-fixed accelerations in terms of the 

rotating coordinate system 1 , l r , 1^, it is convenient to employ the 

c 

vector differential operator on R: 
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FIGURE 1. RE-ENTRY COORDINATE SYSTEMS 
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= ~r~ R + ft X R 
at 


( 1 ) 


where the first term is the derivative of the geocentric displacement 
relative to the space-fixed coordinate system l x , ly, l z> the second term 
is the derivative of the geocentric displacement relative to the rotating 
coordinate system 1 , 1^, 1^, and the third term is the angular velocity 

c 

vector of the rotating coordinate system relative to the space-fixed 
coordinate system. 

Referring to Figure 1, it is apparent that the total angular 
velocity vector may be resolved into two component vectors representing 
angular velocity in the equatorial plane (A) and angular velocity in the 
polar plane (<t> c ) : 

ft = A lz + <t> (- 1.) (2) 

C A 


Resolving the vector representing the angular velocity in the equatorial 
plane still further in the polar plane: 

= A cos 4> 1 + A sin 4> 1 - <t> 1 A (3) 

c <$> c r c A 

c 

Substituting equation 3 into equation 1 and expressing R as rl^ yield 
the familiar space-fixed velocity: 


d 

dt 


R = (rl ) + (A cos O 1, + A sin 4> 1 - 0 1 A ) 

dt r c4> crcA 

s c 


X (rl ) = r 1 + r A cos 0 1. + r <J> 1 

r r c A c 4> 


Thus, the space-fixed component velocities are: 


v = r^ 

0 c 

c 


v, = r A cos <t> 

A c 

A substitution of equation 4 for R in equation 1 yields the desired 
space-fixed accelerations. 


(4) 

(5) 

( 6 ) 
(7) 


5 



dtj dt 


R = — ) 4r ) R + nx- 


dt 


dt 


= TT ) (r <L K + r 1_ + r A cos <t> c 1 A > 


c <t> 


+ (A cos $ 1 A + A sin <t> 1 - i 1.) 

c <p c r c a 

c 

X (r 1. + r 1 + r A cos <t> 1.) 

c 0 r c A 


= r 0 + r 0 1. + r 1 + r A cos 4> 1 A 

c 0 c 0 r c A 

c c 

+ r A cos <t> 1, - r ♦ A sin 0 l.+rA cos 4> 1 A 

c A c c a c 


- r A 2 cos 2 0 1 - r H sin 4> 1 A 

c r c c A 

+ r A 2 sin 4> cos 4> - r 4> 2 1 + r 0 1 


c r 


c 0 


* • o \ 

= (r 0 + 2r 4) + r A 2 sin 4) cos ♦ ) 1 . 

x c c c c 4> 


+ (r - r 0 2 - r A 2 cos 2 4> ) 1 
c c r 


+ (r A cos 4> c + 2r A cos 


0 - 2r 0 A sin 


» ) 

c A 


( 8 ) 


Thus, the space-fixed component accelerations are: 

a =r0 +2r0 + r A 2 sin 0 cos 0 

0 c c c c 

c 


a =r-r0 2 -rA 2 cos 2 0 
r c c 

•• m m • • 

a, = r A cos $ + 2r A cos 4> - 2r $ A sin 4> 

A c c c c 


(9) 

( 10 ) 

( 11 ) 


It is convenient to make the substitution of earth longitude for 
space-fixed longitude so that the equations of motion may then be expressed 
in terms of geocentric latitude, earth longitude, and geocentric altitude, 
which are the usual terrestrial navigational parameters. This is 
accomplished by the angular velocity substitution: 
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A = X + Ol) (12) 

where X is the earth longitude and to is the angular_velocity of the earth, 
(Note that the new rotating coordinate unit vector 1^ is collinear with 

the old rotating coordinate unit vector 1^.) The space-fixed velocity 

and space-fixed acceleration then become: 


— R = r 4> 1 + r 1 + r (X + w) cos 4) 1, 

dty c 0 r c X 

' s c 


— p ) R = [r cj> + 2r 4> + r (A, + go) sin 45 cos 4) ] 1 

dt / c c c c 45 


+ [r - r 4> c 2 - r (A, + to) cos 2 4> c ] 1^ 

• • • • 

+ [r A, cos 0 + 2r (A, + go) cos 45 

c c 

- 2r 4> (X + co) sin 4> ] 1^ (14) 

c c a 

This acceleration equation is now of the form which can be substituted 
into Newton's Second Law in order to obtain the equations of motion for 
a re-entry trajectory. 


SECTION IV. DERIVATION OF GRAVITATIONAL FORCES 


The gravitational potential external to the surface of an oblate 
earth is assumed (Ref. 1) to be: 


TT GM , J ( E \ ( 1 o • 2 A N 

u = - — 1 + 3 ( — ) (1-3 sin 2 4) c ) 


where G is the universal gravitational constant, M is the mass of the 
earth, r is the geocentric displacement of the point in consideration, 

J is the second order figure constant for the earth, r^ is the equatorial 

radius of the earth, and 4> c is the geocentric latitude. 

The gradient vector operator in spherical coordinates is: 


y __ i 1 + — l + i — l 

v d<\> 45 bv r r cos 4> bX X 

c c c 
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The gravitational force per unit mass may be obtained by taking the 
negative of the gradient of the gravitational potential: 



m 


- VU 


i 

r d<l> 


GM 

r 




(1-3 sin 2 $ ) 
c 



3 sin 2 


V-l 




3 \Ty (1 - 3 si ” 2 


sin 2<J> 1 

c <t> 


GM 

„2 




(1-3 sin 2 ♦ ) 
c 


1 * 01 . 
r X 


(17) 


Thus, the gravitational components of acceleration for the earth are: 


G 

♦ 

c 



sin 2<t> 

c 


G =- 

r 


GM 

r 2 



(1 


- 3 sin 2 




(18) 

(19) 

( 20 ) 


These components are the force per unit mass terms caused by the gravita- 
tional acceleration. These terms are to be used in Newton's Second Law. 


SECTION V. DERIVATION OF AERODYNAMIC FORCES 


For the derivation of the aerodynamic forces on the re-entry vehicle, 
it is necessary to depart from the particle dynamics approach long enough 
to adopt a rigid body model of the vehicle by which aerodynamic parameters 
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may be defined. The vehicle is symbolically represented by a blunted 
conical frustum having a sector of a sphere as its blunt base. The 
vehicle model has an auxiliary, controlled, moveable lateral surface by 
which an angle of attack with respect to the vehicle symmetry axis may 
be generated in the vehicle symmetry plane. The vehicle rigid body 
dynamics are assumed not to be coupled with the exterior ballistics of 
the vehicle center of mass, however, and the aerodynamic lift and drag 
forces are considered as acting through the vehicle center of mass 
exclusively. Then, the problem remains essentially one of particle 
dynamics, with the concept of a rigid body vehicle introduced only for 
the purpose of defining the direction of the lift and drag forces with 
respect to the aerodynamic velocity vector. 

The aerodynamic velocity is the velocity of the vehicle center of 
mass with respect to the ambient atmosphere. The earth's atmosphere is 
assumed to be stationary with respect to the rotating earth, with any 
variations being taken into consideration as wind effects. The 
aerodynamic velocity is computed as the relative velocity of the vehicle 
center of mass to a coincident air molecule point mass. The space-fixed 
velocity of the air particle in the stationary atmosphere is assumed to 
be : 


V =0 1 +01 + rco cos 4) 1. (21) 

A„ 4> r c A, 

Cj C 

Assuming that the wind effects of the air mass are restricted to a mass 
flow which is laminar with respect to geocentric altitude, the velocity 
of the air particle then becomes: 

= V w c ° s \\ + o T. + (r w cos $ c + V y sin (22) 

W c 

where and A^. are the speed and azimuth, respectively, of the air flow 

with respect to stationary air at a geocentric displacement r. The 
aerodynamic velocity is then obtained by taking the difference of 
equation 22 from equation 13 to yield the relative velocity: 

? A = <r *c - V W cos V T * + * T r 

c 

+ (r X cos sin A^) 1^ (23) 

For the resolution of the lift and drag vectors into the rotating 
coordinate system, it is convenient to define two reference angles to 
describe the orientation of the aerodynamic velocity vector in that 
coordinate system. Thus, the dive angle of the aerodynamic velocity 
vector is defined to be: 
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-l 

cos 



(24) 


and the azimuth of the aerodynamic velocity vector is defined to be: 


A z = tan 


•i Jli_ 

r 


COS q> 

c 

l - v, 


v w sin \ 


w 


COS 




(25) 


The aerodynamic drag vector is defined to be collinear with the 
aerodynamic velocity vector but oppositely directed. The aerodynamic 
lift vector is defined to be normal to the aerodynamic velocity vector, 
oriented in the vehicle symmetry plane, and directed oppositely from the 
vehicle's auxiliary lateral surface. 


The resolution of the aerodynamic lift and drag vectors into the 
rotating coordinate axes is illustrated in Figure 2. The bottom view 
represents the local geocentric vertical plane which contains the 
aerodynamic velocity vector, and the top view represents the local 
geocentric horizontal plane. It should be clarified here that the 
bank angle y is defined to be the angle between the lift vector and the 
vertical geocentric plane containing the aerodynamic velocity vector. 

The lift and drag resolutions are respectively: 

= L (- sin 7 sin - cos y cos 0 cos A^) 1^ 

u c 

+ L (cos 7 sin © A ) 1 

+ L (sin 7 cos A z - cos y cos 0 sin A z ) 1^ 

F = D (- sin 0. cos A„) 1 + D (- cos 0.) 1 

D A Z q> A r 

c 

+ D (- sin 0^ sin A z ) 1^ 

With the lift and drag directions resolved into the rotating 
coordinate system axes, there remains the magnitudes of those aerodynamic 
forces which are respectively: 

L = q S C L = p V A 2 S (a L a + b L a 2 ) (28) 


(26) 


(27) 


D = q S C D = \ p V/ S (C D + a D a + b D a 2 ) 

o 


( 29 ) 
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where p is the atmospheric mass density, S is the vehicle total effective 
aerodynamic surface area, and a^> b^, , a^, and are empirical 

o 

coefficients which determine the total lift and drag coefficients C.^ 
and (2^, respectively, as a function of aerodynamic angle of attack and 
Mach number. 

Equations 28 and 29 yield the aerodynamic forces to be used in 

: equations of motion. The aerodynamic 

- cos 7 cos © A cos A z ](a L a+b L a 2 ) 

a D a + b D 1 <t> 

c 

(a L a + b L a 2 ) 

b D a 2 )) T. 

+ — {[sin 7 cos A - cos 7 cos 0 A sin A„](a T a+b 7 a 2 ) 

2 m Z A Z L L 

+ [- sin 0 A sin k^} (C D + a D a + b D a 2 )) 1^ (30) 

o 

SECTION VI. DERIVATION OF THE COORDINATE TRANSFORMATION EQUATIONS 


Newton's Second Law for the vehicle 
force per vehicle mass then becomes 


f a pV s ,, . . 4 

r=2~ t[ - sin 7 sin V 


+ [- sin Q ^ cos A z ] (C^ + 

o 


p V s 

+ 2 t [cos 7 sin 0 A ] 


+ [- cos 0 A ] (C D + a D a + 


During the course of the re-entry trajectory, the orientation of the 
rotating coordinate axis system in which the equations of motion are 
expressed is continuously changing with respect to fixed space. It is 
desired to express certain inertial and space-fixed vector quantities of 
the moving coordinate system in terms of the inertial platform coordinate 
system defined in Section II. This resolution may be accomplished by an 
appropriate sequence of Euler angle rotations of the moving coordinate 
axes and a subsequent coordinate translation of the origin of the 
rotated coordinate system from the center of the earth to the surface of 
the earth. 
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This Euler angle sequence is illustrated in Figure 3. The first 
Euler rotation is about the instantaneous 1^ axis through an angle 

equal to the instantaneous geocentric latitude. The resolution equations 
expressed in matrix form are: 


1 


cos <t> sin 4> 0 


i 4 

X 1 


c c 


*0 

I 


- sin <t> cos <t> 0 


1 



c c 


r 

I 

_ z l_ 


0 0 1 


h 


The second Euler rotation is about the newly generated lx^ axis through 
an angle equal to the space-fixed longitude. The resolution equations 
expressed in matrix form are: 



1 


y z 


1 


1 0 

0 cos A 

0 sin A 


sin A 
cos A 


(32) 


The third Euler rotation is about the newly generated lz 2 axis through 
an angle equal to the initial geocentric latitude. The resolution 
equations expressed in matrix form are: 




o 


- 



cos 4> - sin <t> 0 


1 

c c 


X 

o o 


2 

sin 4> cos 4> 0 


1 

c c 

o o 


y 2 

0 0 1 


1 



Z 



2 


(33) 


The fourth Euler rotation is about the initial l r axis through an angle 
equal to the initial space-fixed azimuth. The resolution equations 
expressed in matrix form are: 
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X 


cos Ap 0 sin A ^ 

1 


0 10 

n 



\ 


- sin Aj 0 cos Aj 



_ — 


( 34 ) 


The total vector resolution is then obtained by performing the appropriate 
matrix multiplications which yield; 


cosA t cos4> cos4> + cosA T cos<l> sin<t> - 

lcc Icc 

o o 

cosA-rSinO cosAsin<t> cosA T sin<t> cosAcos<t> cosA T sin<t> sinA 

Ic cic c ic 


-sinA T sinAsin<l> 
l c 


sm 4> cos4> 
c c 
o 

-cos<f>c cosAsin<t> 
o c 


+sinA- r sinAcos<!> 

I c 


sm4> s m<t> + 
c c 
o 

cost J> cosAcos<t> 
c c 

o 


H-sinA^cosA 


-cos<t> sinA 
c 

o 


-sinA T cos4> cos<J> + -sinA T cos<t> sin4> + 


sinA-rSin^ cosAsinO sinA T sin<l> cosAcos<I> -sinA T sin<i> sinA 
Ic cic c ic 


-cosA T sinAs inty 
I c 


+cosA,sinAcos4> 

I c 


+cosA^.cosA 


( 35 ) 


This combined matrix will then transform the space-fixed displacement, 
velocity, or acceleration vectors from the instantaneous space-fixed 
orientation of the moving coordinate system into the idealized inertial 
platform coordinate system orientation. A subsequent translation of the 
origin along the 1 axis then yields the final idealized inertial 

platform coordinate system, i.e., 


rip = n E 


r 

o 


( 36 ) 


where 7]p is the platform rj, is the resolved space-fixed r), and r Q 
is the radius of the earth at the re-entry point. 
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SECTION VII. ASSEMBLAGE OF THE FINAL EQUATIONS 
OF MOTION FOR NUMERICAL SOLUTION 


The equations of motion may now be obtained by equating the results 
of equation 8 to the results of equation 17 plus equation 30. Cancellation 
of the unit vectors then yields the three simultaneous scalar differential 
equations : 

2 


r d> + 2r <t> + r (X + w) sin <t> cos <t> 

c c c c 




sin 24> 


- ” sin 7 sin A„ - — cos y cos 9. cos A, - - sin 0. cos A„ (37) 

IQ « m A £ m A L 


" is /■>* \2 o , GM 

r - r 4> - r (X + w) cos <t> = - ~ 


(1-3 sin 2 <t> c ) 


L , _ D _ 

+ “ cos 7 sm - — cos B A 
m Am A 
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r X cos <t> + 2r (X + w) cos <t> - 2r <t> (X + co) sin <)> 

c c c c 


= — sin 7 cos A„ - — cos y cos 9. sin A„ 
m ' Z m A Z 


- — sin 9. sin A„ 
m A Z 


(39) 


Solving each of the differential equations for the highest derivative 
then puts them in a form solvable by numerical integration on the IBM 7090 
digital computer: 
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— - (X + u) 2 sin <t> cos <t> - J sin 2<P 


c r \r / 
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r = r 0 2 + r (X + w) cos 2 0 
c c 


GM 

v.2 


1 + J ( — ) (1-3 sin 2 0 ) 

v r J c 


+ — cos 7 sin 9 . 
m ' A 


- — cos 9 
m A 


X = - — — + ^ +20 (X + go) tan 0 
r c c 

1 fh . L „ . A 

+ — -S“ sin 7 cos A„ - — cos 7 cos sm A 

r cos 0 m Z m A Z 


( 41 ) 


- - sin 9 sin A \ (42) 

m A Z J 

To solve this set of differential equations, it is necessary to 
provide six initial conditions. Three of those initial conditions may 
be vspecif ied by considering the initial position of the vehicle. The 
initial geocentric latitude and initial earth longitude may be specified 
directly. The displacement may be derived by specifying the initial 
geocentric altitude and applying the formula derived in Figure 4. 



(43) 


The remaining three initial conditions may be derived by considerations of 
the initial space-fixed velocity of the vehicle: 


r 0 = V_ sm 9- cos A_ 

o c I I I 

o 


r = V_ cos 9 T 
o I I 


r (X + go) cos 0 = V T sin 9 r sm A T 

00 cl I I 

o 


(44) 

(45) 

(46) 


where V^. is the initial space-fixed velocity, 9 ^ is the initial space- 

fixed dive angle, and A^. is the initial space-fixed azimuth. The initial 

space-fixed azimuth might also be expressed in terms of the inclination 
(j) of the vehicle's instantaneous orbit plane to the equatorial plane: 
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The initial- conditions to be obtained from the initial velocity 
considerations thus become: 

V_ sin 6 _ cos A T 

: I I I 


( 47 ) 


(48) 


r 

o 



cos Oj 


X 

o 


Vj sin sin 


r cos <t> 
o c 

o 


(49) 

(50) 


SECTION VIII. INPUT REQUIREMENTS FOR THE COMPUTATIONAL PROCEDURE 


The input requirements for the £BM 7090 re-entry program fall broadly 
into three categories: (a) geophysical data (b) aerodynamic data, and 
(c) initial conditions for the equations of motion. 


The geophysical data required are the Gaussian constant (GM) , second 
order figure constant (J), the equatorial radius (r E ), the polar radius 


(r^) , and the angular velocity of rotation (w) . 


The first four parameters 


may be manipulated to produce any oblate spheroid and corresponding 
gravitational field desired. The values used by the author for those 
five items can be found in Reference 2. 


The aerodynamic data required are atmospheric phenomena and vehicle 
characteristics. The Computation Division provides 1959 ARDC standard 
atmospheric mass density and speed of sound data, both as a function of 
altitude. Provision is made to generate any variation from a standard 
atmosphere by a density modification factor which is constructed as a 
curve read function of altitude. A wind velocity profile as a function 
of altitude and a constant wind azimuth are also available. Provision 
is made for a constant, effective aerodynamic surface area to vehicle mass 
rates (S/m), a constant angle of attack (a), and a constant bank angle 
( 7 ). The empirical constants (a^, b^, c^ , a^, and b^) which determine 

o 

the lift and drag coefficients may be specified as a curve r^ad function 
of Mach number. 
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The initial conditions of the equations of motion may be effected 
via Section VI by specifying the vehicle's initial space-fixed velocity 
(Vj), initial space-fixed dive angle (0^.), initial space-fixed azimuth 

(Aj) or instantaneous orbital inclination (j), initial earth longitude 

(X ) , initial geocentric latitude (4> ) , and initial geocentric altitude 

o 

<*>„>• 


The print-out time interval (At) may be chosen to be any multiple 
of the integration step time interval 1 ^ (At) which is also arbitrary accord- 
ing to the computational accuracy desired. 

SECTION IX. CONCLUSION 

A system of equations has been derived which can be solved numerically 
on the IBM 7090 digital computer to describe the motion of a vehicle re- 
entering the atmosphere of the earth. Components of acceleration, velocity, 
and displacement as a function of time can be obtained in a spare-fixed 
coordinate system which is parallel to the space-fixed platform on board 
the vehicle. Thus, the trajectory simulation given by these equations is 
easily adaptable to hardware accuracy studies. 
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